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Abstract 

We compute the matrix elements of the operators which contribute 
to spectator effects in inclusive decays of i?-mesons. The results agree 
well with estimates based on the vacuum saturation (factorization) 
hypothesis. For the ratio of lifetimes of charged and neutral mesons 
we find T{B~)/T{Bd) = 1.03 it 0.02 it 0.03, where the first error rep- 
resents the uncertainty in our evaluation of the matrix elements, and 
the second is an estimate of the uncertainty due to the fact that the 
Wilson coefficient functions have only been evaluated at tree-level in 
perturbation theory. This result is in agreement with the experimen- 
tal measurement. We also discuss the implications of our results for 
the semileptonic branching ratio and the charm yield. 



1 Introduction 

Inclusive decays of heavy hadrons can be studied in the framework of the 
heavy quark expansion, in which widths and lifetimes are computed as series 
in inverse powers of the mass of the 6-quark H, ^ (for recent reviews and 
additional references see refs. @, ^). The leading term of this expansion cor- 
responds to the decay of a free-quark and is universal, contributing equally 
to the lifetimes of all beauty hadrons. Remarkably there are no corrections of 
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0(l/mfe), and the first corrections are of 0(1 /ml) ^. "Spectator effects", 
i.e. contributions from decays in which a hght constituent quark also partici- 
pates in the weak process, enter at third order in the heavy-quark expansion, 
i.e. at 0{l/ml). However, as a result of the enhancement of the phase space 
for 2^2 body reactions, relative to 1— >3 body decays, the spectator effects 
are likely to be larger than estimates based purely on power counting, and 
may well be significant. The need to evaluate the spectator effects is rein- 
forced by the striking discrepancy between the experimental result for the 
ratio of lifetimes M : 



r(A, 



b 



r(5, 

and the theoretical prediction 



0.78 ±0.04, (1) 



0.98 + O(l/m^) . (2) 



In order to explain this discrepancy in the conventional approach, the higher 
order terms in the heavy quark expansion, and the spectator effects in par- 
ticular, would have to be surprisingly large. In this paper we compute the 
matrix elements of the four-quark operators needed to evaluate these spec- 
tator effects for mesons; a computation of the effects for the A;, is in progress 
and the results will be reported in a future publication. 

The experimental value of the ratio of lifetimes of the charged and neutral 
S-mesons is [0 

^ = 1.06 ±0.04, (3) 

to be compared to the theoretical prediction 

t(B-) 



l + 0{l/mt). (4) 



Below we determine the contribution to the 0(l/m^) term on the right-hand 
side of eq.(^ coming from spectator effects, which we believe to be the largest 
component. The same spectator effects also contribute to the semileptonic 
branching ratio of the i?-mesons, and to the charm-yield {ric, the average 
number of charmed particles in decays of 5- mesons). This will be briefly 
discussed in section |. 
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At 0(1 /ml), the non-perturbative contributions to the spectator effects 
are contained in the matrix elements of the four-quark operators: 

O'y.A = h-f,qLqLrbL (5) 

Ol_p = hRqiqihR (6) 

= hl.T'^qLqiYT^L (7) 

T|_p = hRT'^qLClLT'-hR , (8) 

where q represents the field of the light quark, and the T"'s are the generators 
of the colour group. Throughout this paper we take these operators to be 
defined at a renormalization scale m;,. The subscripts L and R represent 
"left" and "right" respectively. For mesons, following ref. 0, we introduce 
the parameters B12 and £1,2: 

^ {B,\0\,_^\B,) ^ ^^V^Si, (9) 



2mB, 



{B,\Ol_p\B,) ^ '-^^B,, (10) 



^ {B,\n_A\B,) ^ ^^i^ei, (11) 



2mB, 



^ 'B,\TI_,\B,) ^ f^^e,, (12) 



2mH 
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where /s^ is the leptonic decay constant of the meson Bq (in a normalization 
in which f^^ ~ 131 MeV). The introduction of the parameters -81,2 and £1,2 is 
motivated by the vacuum saturation (or factorization) approximation |jlO|| in 
which the matrix elements of four-quark operators are evaluated by inserting 
the vacuum inside the current products. This leads to i?j = 1 and = 
at some renormalization scale, /x, for the operators. It may be argued that 
the scale at which the factorization approximation holds is different from 
our chosen renormalization scale m^, specifically that if the approximation is 



valid at all then it should hold at a typical hadronic scale of about 1 GeV |TT 
In QCD, in the limit of a large number of colours 

Bi = 0(1) and Ei = 0{1/N,) . (13) 

In this paper we evaluate the parameters B^ and Ei in lattice simulations, 
thus testing the validity of the factorization hypothesis. 
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The results of our calculations indicate that the vacuum saturation hy- 
pothesis is (surprisingly?) well satisfied. In the MS renormalization scheme 
we find: 



Bi{mb) = 1.06(8) B2imb) = 1.01(6) (14) 
£i(mb) = -0.01(3) e2{mb) = -0.01(2) (15) 



and 

t(B' 



1.03 ± 0.02 ± 0.03 , (16) 



in good agreement with the experimental value in eq. (^. 

The present calculation is very similar to that of the i^^-parameter of B-B 
mixing, for which several recent simulations have been performed |jT2|, [13|, [14 



including one using the same configurations used in this study [0. We 
use the calculations of the Bb parameter as a comparison and check on our 
results, both for the perturbative matching coefficients and for the evaluation 
of the matrix elements. In performing this comparison we believe that we 
have found an error in the perturbative calculation. We also stress that the 
feature that the values of the matrix elements of the operators (^-(^ are 
close to those expected from the vacuum saturation hypothesis is also present 
in the evaluation of the iJ^-paraiiieter. 

The plan of the remainder of this paper is as follows. In the next section 
we derive the relation, at one-loop order of perturbation theory, between the 
operators defined in eqs. (^-(||) and the bare lattice operators whose matrix 
elements are computed directly in lattice simulations. The technical details 
of the calculation are relegated to appendix A. In section ^ we present a 
description of the lattice computation, the results of this computation and 
a discussion of the implications. The comparison of our study to that of 
B^-B^ mixing is presented in section ^ (the discussion of the evaluation of 
the matching coefficients for B-B mixing is described in detail in appendix 
B). Finally, in section |]we present our conclusions. 

2 Perturbative Matching 

The Wilson coefficient functions of the operators (^-(i) in the OPE for in- 
clusive decay rates have been evaluated only at tree- level Pj. At this level of 
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precision it is sufficient to compute the matrix elements in any "reasonable" 
renormalization scheme. In this paper we present the results for the matrix 
elements of the HQET operators in the MS* scheme at a renormalization 
scale rrih [|. In order to obtain these from the matrix elements of bare oper- 
ators in the lattice theory, with cut-off a^^ (where a is the lattice spacing), 
which we compute directly in our simulations, we require the corresponding 
matching coefficients. In this section we present these matching coefficients 
(at one-loop order in perturbation theory); we postpone a detailed descrip- 
tion of their evaluation to Appendix A. In the following section we will use 
these coefficients and the computed values of the matrix elements of the bare 
operators to determine the coefficients Bi and in the M^* scheme. 

The coefficients presented in the section were all obtained using local 



lattice operators defined in the tree-level improved SW-action |jT5[ defined in 
eq. (piD below. 



It is convenient to perform the matching in two steps: 

i) The ffist step is the evaluation of the coefficients which relate the HQET 
operators in the continuum {MS) and lattice schemes, both defined at 
the scale a~^, 

0?{a-') = Of {a-') + ^a.,Of{a-') , (17) 

where the lattice (continuum) operators are labelled by the superfix L 
(C). The mixing of operators is such that it is necessary to compute the 
matrix elements of eight lattice operators Oj'(a~^) (see table |1|, where 
the coefficients Dij are presented). 

ii) We then evolve the HQET operators in the continuum scheme from 
renormalization scale fi = a^^ to scale fi = m^. 

Of (mt) = M,,{a-\mi,)Ofia-^) . (18) 

This evolution, sometimes known as "hybrid renormalization" , requires 
knowledge of the anomalous dimension matrix [|l^, 

The matching procedure involves short-distance physics only, and so can 
be carried out in perturbation theory. Lattice perturbation theory generally 



"'^From these it is possible to obtain the matrix elements in any other renormalization 
scheme using perturbation theory. 
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Dij 


i = 1 


i = 2 


i = 3 


i = 4 


O 


J = l 


-21.64 


- 


2.06 


0.54 




J = 2 




-21.64 


2.16 


2.40 




J = 3 


9.29 


2.43 


-10.80 


2.83 




J = 4 


9.72 


10.79 


11.34 


-9.05 




J = 5 


-18.37 




-3.06 




bRYQR qii^h 


j = Q 


36.75 


18.37 


6.12 


3.06 




J = 7 


-13.78 




6.89 






J = 8 


27.56 


13.78 


-13.78 


-6.89 


IrT-Ql qRT-bL 



Table 1: Coefficients Dij for the four operators of eqs.(H)-(||). 



converges very slowly, and therefore, where possible, the matching should be 
performed non-perturbatively so as to minimise these systematic errors. In 
this letter, however, we do not use a non-perturbative renormalization, but, 
when evaluating the matching coefficients, we do use a "boosted" lattice cou- 
pling constant in order to partially resum the large higher order contributions 
(e.g. those coming from tadpole graphs). 

In order to obtain the parameters Bi and Ei it is also necessary to deter- 



mine the normalization of the axial current [1^ . In this case there is a single 
coefficient defined by 

{0\A^{a-')\B)=Z^{0\A^{a~')\B) , (19) 

which at one loop order in perturbation theory is given by: 

Za = 1 + (4 - - X, + X3 j ^ 1 - 20.0 = 0.79 (20) 

where the coefficients Xj are defined in eq. (^) and tabulated in table |^ in 
Appendix A. Aq and A^ are the time components on the axial current in the 
HQET in the MS" and lattice schemes respectively, both defined at the scale 
a~^. Using hybrid renormalization one can obtain the axial current in the 
MS* scheme at scale mf,. We stress again that the results presented for the 
-Bj's and ej's below were obtained with both the four quark operators and 
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the axial current defined in the HQET in the MS" scheme at the scale rrih 0. 

In the following section we combine the results of the matrix elements 
computed on the lattice with the perturbative coefficients presented in this 
section to obtain the -Bj's and e^'s. 



3 Lattice Computation and Results 

The non-perturbative strong interaction effects in spectator contributions to 
inclusive decays are contained in the matrix elements of the eight four-quark 
operators, Oj, given in table |l]. We evaluate these matrix elements in a 
quenched simulation on a 24^ x 48 lattice at /3 = 6.2 using the SW tree-level 



improved action |T5 



gsw ^ ^gauge ^ ^Wilson - ^ q{x) F,,{x)a q{x) , (21) 

where 5'^'^"^'^ and 5'^^^'^°'^ are the Wilson gauge and quark actions respectively. 
The use of this action reduces the errors due to the granularity of the lat- 
tice to ones of O^aso), where a is the lattice spacing. We use the 60 SU(3) 
gauge-field configurations, and the light quark propagators corresponding to 
hopping parameters k = 0.14144, 0.14226 and 0.14262 which have been used 
previously to obtain the 5-parameter of B^-B^ mixing |T^ and other quan- 
tities required for studies of 5-physics. The value of the hopping parameter 
in the chiral limit is given by Kc = 0.14315. The calculation of the matrix 
elements of the operators in table |^ is very similar to that of the AB = 2 
operators from which the 5-parameter is extracted, and we exploit the sim- 
ilarities to perform checks on our calculations (see section ^ and Appendix 
B). 

The evaluation of the matrix elements requires the computation of two- 
and three-point correlation functions: 

C(t.) = 5^(0|J(x)jt(o)|0) , (22) 

where we have assumed that > 0, and 

K{t,, ty)^Yl{0 I J{y)Of{0)j\x) I ) , (23) 

x,y 



^These differ at one-loop level from the corresponding operators in QCD. 
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where ty > > t^. In eqs. P^ ) and ( P^ and J are interpolating operators 
which can create or destroy a heavy pseudoscalar meson (containing a static 
heavy quark) and Oj is one of the operators whose matrix element we wish 
to evaluate. In practice we choose J to be the fourth component of an axial 
current. It is generally advantageous to "smear" the interpolating operators 
J and over the spatial coordinates, in order to enhance the overlap with the 
ground state. Following ref. |jl2| we have used several methods of smearing 



and checked that the results for the matrix elements of Oj are independent of 
the choice of smearing. In this paper we present as our "best" results those 
obtained using a gauge-invariant smearing based on the Jacobi algorithm 
described in ref. [20|. The smeared heavy-quark field at time t, b^{x,t), is 
defined by 



^M{x,x')b{x',t) 



(24) 



where 



N 

M{x,x') = ^KsA"(f,f') 

ra=0 



(25) 



and A is the three-dimensional version of the Laplace operator. The param- 
eter Ks and the number of iterations can be used to control the smearing 
radius, and here we present our results for = 0.25 and = 140 which 
corresponds to an rms smearing radius tq = 6.4a |2l]|. The smeared interpo- 
lating operator is then chosen to be J^{x) = W {x)'~^^q{x) , where q is the field 
of the light quark. The local interpolating operator is simply defined to be 
J^(x) = 6(x)75g(x). 

The eight matrix elements {B\ Oj^(O) \B) are determined by computing 
the ratios: 



Rj(ti,t2 



(26) 



where ti and t2 are positive and sufficiently large to ensure that only the 
ground state contributes to the correlation functions. The indices S and 
L denote whether the interpolating operators are smeared or local. It is 
convenient to choose both J and to be smeared in the three-point function 
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Kj and to evaluate the two-point functions with a local operator at the source 
and a smeared one at the sink. At large time separations ti and t2, 



Rj{tuh)^-—^{B\0!fiO)\B) , (27) 



where Zl is given by 



(0| J^(0) \B) = V2rnBZL . (28) 



The Rj defined in eg. (|26|) contribute directly to the _Bj's and the e^'s, 



apart from perturbative matching factors. To see this, note that the leptonic 
decay constant of the 5-meson in the static limit (i.e. infinite mass limit for 
the fe-quark) is given by 

/I = , (29) 



so that 



TUB 



AZ2 

R,{h, h) - -^{B\ Of{0) \B) = -j^{B\ Of{0) \B) , (30) 

^B^L Jb^B 

which corresponds precisely to the normalization of the Sj's and si's in 
eqs. (PI)-(P^) (apart from the matching factors Z\ and that of the four-quark 
operator Oj described in section ^. 

In these computations it is particularly important to establish that the 
contribution from the ground-state has been isolated. The natural way to do 
this is to look for plateaus, i.e. for regions in ti and ^2 for which Rj{ti, ^2) is 
independent of ti and ^2- Since the statistical errors grow fairly quickly with 
ti^2 our ability to verify the existence of plateaus is limited. For example in 
fig. |I] we present our results for Rj(ti, ^2), obtained with the light quark mass 
corresponding to k = 0.14226 as a function of ^2 for ti = 3. The results are 
consistent with being constant, but the errors are uncomfortably large for 
^2 > 5 or so. An alternative, and perhaps more convincing, way to confirm 
that the contribution from the ground-state has been isolated is to check 
that the values of the Rj are independent of the method of smearing used to 
define the smeared interpolating operators. Following ref. [^, in addition 
to the gauge invariant prescription for smearing described above, we have 
defined the smeared field in four other ways (having transformed the fields 
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Lrauge inv. 


Exponential 


vjraussian 


Uube 


Double Cube 


i?i(3,3) 


1.02(3) 


1.03(2) 


1.04(3) 


1.03(2) 


1.07(8) 


i?2(3,3) 


1.00(2) 


1.01(2) 


1.01(3) 


1.01(2) 


1.00(5) 


i?3(3,3) 


0.00(1) 


-0.00(1) 


-0.00(2) 


-0.00(1) 


-0.03(3) 


-R4(3, 3) 


0.00(1) 


0.00(1) 


0.01(1) 


0.00(1) 


0.01(2) 


i?5(3, 3) 


-0.99(4) 


-1.01(3) 


-1.01(4) 


-1.00(3) 


-1.01(7) 


-Re (3, 3) 


-0.99(3) 


-1.00(2) 


-1.01(3) 


-1.00(3) 


-1.01(6) 


^7(3,3) 


-0.04(2) 


-0.03(2) 


-0.04(2) 


-0.03(2) 


-0.09(5) 


-^8(3, 3) 


-0.00(1) 


-0.00(1) 


-0.01(1) 


-0.01(1) 


-0.02(2) 



Table 2: Values of the Rj{3, 3) obtained with different smearing methods for 
a light quark with k = 0.14226. 



to the Coulomb gauge). For these additional four cases, M{x,x') of eq. (pSQ 
is replaced by the following: 

Exponential : M{x,x') = exp(— |x — x'\/ro), (31) 

Gaussian: M{x,x') = exp{—\x — x'\'^/rl), (32) 

3 

Cube: M{x,x')=Y[e{ro-\xi-x[\), (33) 

1=1 

Double Cube: M(x,x')= TT - ~ ^'^ ^ Q(2ro - \xi - x'A), (34) 

and we have chosen = 5. As an example, we present in table the results 
for the Rj{ti, for ti = 3 and t2 = 3, again for the middle value of the three 
ft's, K = 0.14226. 

We take as our best results those obtained with the gauge invariant smear- 
ing, and with ^1 = ^2 = 3. After extrapolating to the chiral limit we find the 
following results for the Rj: 

i?i(3,3) = 1.04±0.04 i?2(3,3) = 1.00±0.03 (35) 

i?3(3,3) = -0.01 ±0.02 i?4(3,3) = -O.OOiO.Ol (36) 

i?5(3,3) = -0.97±0.06 i?6(3,3) = -0.98±0.05 (37) 

i?7(3,3) = -0.03±0.03 i?8(3,3) = -0.01 ±0.01 (38) 

We now combine these results with the matching coefficients presented 
in section to determine the i?,'s and the £,'s. For the value of the lattice 



10 




258 11 258 11 258 11 258 11 



hLYT'^qmLl^T^L hRT'^qL-qiJ^R hR-i^'T''qRqi^-ii'T%L bRT^qLqRT%L 
0.5 I — I I I I I I I I I I — I I — I I I I I I I I I I — I I — I I I I I I I I I 



< 0.0 

CO 



-0.5 



••••• 







2 5 8 11 




2 5 8 11 2 5 8 11 
distance ^2 (^i = 3) 



I I I I I I I I I 











2 5 8 11 



Figure 1: Results for the j = 1 ■ ■ ■ 8 as a function of t2 for ti = 3. The 
value of the quark mass corresponds to k, — 0.14226. 



coupling constant we take one of the standard definitions of the boosted 
coupling: 

—A = TT^ = 0-0105 . 39 
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Using the coefficients from section ^ we then obtain 



{Zj,)-^[l + ^^^D^ i?,= 0.98(8 



B,.{a-^\ 



{ZaY 



{ZaY 



e2{a 



v-2 



1 + 
1 + 
1 + 



as 


(a 


') 




Att 




as 


(a" 


') 




Arc 




as 


{a' 


') 




47r 




as 


(a" 




Atc 



d] Rj = 0.93(5) 
L>) i?j= 0.01(3) 
d] Rj = 0.00(2) 

/ 4j 



We now evolve those coefficients to the scale nib, using the relations: 



1 + 



2Cf6 



B(a 



25 



■eia 



1 + 



eia — 



1^ ^^'^Efa-^: 



iV2 



where 



agio, 
a 



9/2/3o 



0.09(3) . 



(40) 
(41) 

(42) 



In estimating S and its uncertainty we have allowed for a conservative 
variation of the parameters around the "central" values {Aqcd = 250 MeV, 
a-^ = 2.9 GeV, rub = 4.5 GeV and = 9). We finally obtain: 



Si(mfe) = 1.06(8) B^inib) = 1.01(6) 
eiimt) = -0.01(3) e2{mb) = -0.01(2) . 



(43) 
(44) 



These matrix elements have also been evaluated using QCD sum-rules [^2] . 
These authors find Bi{mb) = 0.96(4), B2imb) = 0.95(2), £i{mb) = -0.14(1) 
and e2{Tnb) = —0.08(1), differing somewhat (particularly for the e^'s) from 
our values. 

Using the results in eqs. (^) and ( ^4]) we obtain the following value for 
the ratio of lifetimes for the neutral and charged i?-mesons: 



r{B-) 
r{B,) 



1 + kiBi + k2B2 + hei + = 1-03 ± 0.02 ± 0.03 



(45) 
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where the coefficients ki are taken from ref. 0. The ffist error in eq. 
is from the uncertainty in the values of the matrix elements in eqs. 



and (P^, whereas the second is an estimate of the uncertainty due to our 
ignorance of the one-loop contribution to the Wilson coefficient function in 
the OPE (this was estimated by varying the matching scale from ■mb/2 to 
2mfc using the procedure described in ref. 0). The value in eq. ( ^5]) is in 
good agreement with the experimental results in eq. (|^). 

There has been a considerable amount of discussion lately as to whether 
the theoretical predictions for the semileptonic branching ratio of the B- 
meson (Bsl) and the charm yield are consistent with experimental mea- 



surements 0, 0. Here we limit our discussion to a comment on the 
imphcations of our results on these two physical quantities. The spectator 
contributions to Bsl and ric take the form: 

ABsL,,pec = biBi + b2B2 + hei + he2 (46) 
Anc,spec = niBi + ng^a + n^ei + 71^62 , (47) 

where the coefficients hi and rii at tree-level are presented explicitly in ref. 
(see also refs. [0, |TB|). The uncertainties in the predictions for Bsl and ra^ 
due to the errors in the matrix elements in eqs. (^31) and ( ^l)) are small (about 
0.1% for both Bsl and ric). Larger uncertainties are due to our ignorance of 
the one-loop contribution to the Wilson coefficient function in the OPE. For 
example, following the procedure described in ref. with a matching scale 
11 = mb/2 we find ABsL,spec = 0.3(1)% and Anc^spec = 0.5(2)%, whereas for a 
scale 2mB the corresponding results are ABsL,spec = —0.1(1)% and An^spec = 
0.0(1)%. For the charm yield these corrections are negligible as expected, and 
for the semileptonic branching ratio they are small. Nevertheless, it would 
be useful to know the one-loop contribution to the coefficient functions in 
order to eliminate the variation obtained by changing the matching scale. 



4 B-B Mixing 

In this section we revisit the phenomenologically important process of B- 
B mixing. The computation of the matrix elements of the relevant lattice 
operators on the same gauge configurations has already been presented in 
ref. [^], and we do not add to these lattice computations. We do, however, 
wish to make two observations: 
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i) The matrix elements of the lattice AB = 2 operators relevant for B-B 
mixing factorize with a similar precision to that found for the four- 
quark operators considered in section ^. Because of the way in which 
lattice computations are usually organised, this property is not as read- 
ily manifest for the AB=2 operators as it is for the spectator effects. 
However, by using colour and spin Fierz identities, we demonstrate that 
the values of the matrix elements of the lattice AB = 2 operators are 
very close to those expected using factorization. Some of our observa- 



tions have already been noted in ref . , where the Wilson formulation 
for the light quarks was used. We extend this investigation of factor- 
ization, and show that each contribution to the matrix elements (i.e. 
each Wick contraction) is close to the estimated value obtained using 
the factorization and vacuum saturation hypothesis. 

ii) We believe that there is an error in the published value of the matching 
factors for the AB = 2 operators using the SW action for the light 
quarks. We discuss this in some detail in Appendix B; in this section 
we briefly comment on the consequences of this error. 

We now consider these two observations in turn. 



4.1 Factorization 

In order to determine the 5-parameter of B-B mixing we need to evaluate 
the matrix elements of the three lattice operators Ol, Or and On defined in 
eqs.(BB)-(pS|) of Appendix B, as well as that of 



Os = hi ki ■ (48) 

We now show that these matrix elements (and related ones) are reproduced 
remarkably accurately by assuming the factorization hypothesis and vacuum 
saturation. For each of these operators (Oj) we define the ratio Rj{ti,t2) 
analogously to eq.(^) as follows: 

E(f f)-^ Kfi^^M) 



where the correlation functions C and Kj are defined in eqs. (p2D and (psf), 
except that in the Kj the operators Oj are now AB = 2 operators, and the 
interpolating operators destroy a i?-meson and create a i?-meson. 
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We start by explaining explicitly what we mean by factorization (com- 
bined with vacuum saturation). The evaluation of the B-parameter requires 
the computation of three-point functions Kj{ti,t2), which are of the formQ: 

( \q{yh'b{y) 6(0)rg(0) 6(0)fg(0) q{xh'b{x) | ) . (50) 

There are four Wick contractions which contribute to the correlation func- 
tion, and it is convenient to track these by introducing a fictitious quantum 
number, labelled by an integer suffix on each field, so that the correlation 
function is written in the form 

{0\q2{yh\{y) Op_f(0) q4{xh%{x)\0) , (51) 

where 

Opf = &irg2 ^sf g4 + h^q4 ^if g2 - h'^q^ hfq2 - hTq2 bitq^ , (52) 

and r,r are Dirac matrices. The contraction of spinor {a, (3) and colour (a) 
indices in each bilinear in eq. (^) is implicit (e.g. biTq2 = ^i, a^a,/35'2,/3)• 
Only contractions between fields with the same suffix are allowed, and the 
four terms in eq. ( |52D correspond to the four original Wick contractions. For 
all the operators of interest, by using colour and spin Fierz identities, it is 
possible to rewrite each of the four terms into sums of operators of the form 
(biT'q2) (bsT'q^) and (biT'T"-q2) {b^T'T^q^), for some 7-matrices T' and T'. In 
the factorization and vacuum saturation hypothesis 

(B|(6irr'^g2)(&3f'r'^g4)|5) - , (53) 

and 

( B I {hr'q2) ihr'q,) \B)c^{B\ {hr'q2) | 0) ( | (ht'q^) \B) . (54) 

Lorentz invariance implies that each of the matrix elements on the right-hand 
side of eq.(l5^) vanishes or is proportional to the leptonic decay constant fs- 
We now consider each of the operators in turn: 

• Ol and Or. The factor of 3/8 in eq.(|i9|) was chosen so that the factor- 
ization hypothesis gives Rl = Rr = 1- Numerically we find: 

RLi3,3) + Rr(3,3) , , 

' ^ . ; ^ 0.95 ±0.03 . (55) 



''The extension of this discussion to include smeared interpolating operators is com- 
pletely straightforward. 
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• Os'- Lorentz invariance implies that the matrix element {B \ba'^'^{l — 
7^)g I ) vanishes. Using this fact we deduce that factorization implies 
that Rs ~ 5/8. The numerical result for Rs is 0.60(3), in very good 
agreement with the estimate based on factorization. 

• Ojv: Finally factorization implies that Rj^i ~ 1, in good agreement with 
the numerical value 0.97(4). 

In order to illustrate further that the numerical results quoted above are 
in agreement with expectations based on factorization and vacuum saturation 
we present separately in fig. |^ the contributions to each of the ratios from 
the two independent contractions: 

R- = {B\{hr'q2){ht'q,)\B) (56) 

and 

R'={B\ihrq,)ihsT'q2)\B) , (57) 
as well as the total combination 

R = h2R'' - 2R'') (58) 
8 

We see that not only are the values of the i?'s quoted above in agreement 
with expectations from factorization, but the values of each of the i?" and 
R'' are also those we would expect on this assumption: 

R^L^-^Rl ; R's^-1 R's^^R's (59) 
R'n^-lRl ; R% ^2 R%^Ir%. (60) 

In this subsection we have demonstrated that the surprising precision of 
predictions for matrix elements obtained using the factorization hypothesis 
in inclusive decays, which was discussed in section |], also applies to existing 
results for B-B mixing. 

4.2 Matching 

In this subsection we discuss very briefly the implications of the error in 
the published value of the perturbative matching factors in B-B mixing (see 
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Figure 2: Results for the i?", i?^ and Rj (j = L, 5, i?, A^) as a function of t2 
for ti = 3. The value of the quark mass corresponds to k = 0.14226. 



appendix B). For example, we take method (a) of ref. [|T2|, in which each of 
the ratios Ri is fitted separately, and find 



B^{m,) = 0.66(2) , = af'^B.im^) = 0.96(3) , 



(61) 



where Bb(mb) and Bf, are the 5-parameter in the MS* scheme at scale 
and the renormalization group invariant i?-parameter respectively. For the 
comparison we only quote the statistical error. In order to see the effect of the 
error in the matching coefficient, the result in eq. (|6l|) should be compared 



with that obtained in ref. |]T2| using an identical procedure (except for the 
values of the matching coefficients): 

B,im,) = 0.69(2) , B, = af'^''B,{m,) = 1.02(3) . (62) 

Gimenez and Martinelli had pointed out that the authors of ref. |T^ 
had used matching coefficient which did not include the contributions to the 
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O(a^) terms in the lattice operators which they had used \T^\. As explained 
in ref. [|l^], this leads to a negligible correction to the 5-parameter. The 
differences in the values in eqs. ( [61]) and (^) is therefore largely due to the 
error discussed in appendix B. 



5 Conclusions 

In this paper we have evaluated the matrix elements which contain the non- 
perturbative QCD effects in the spectator contributions to inclusive decays 
of i?-mesons. For the MS scheme our principal results are contained in 
eqs. ( ^31) and (0). The raw lattice results from which the matrix elements 



in any renormalization scheme can be determined are given in eqs. (|35D -(p8D. 

The results for the matrix elements are very close to those which would be 
expected on the basis of factorization and the vacuum saturation hypothesis 
(particularly for the matrix elements of the bare lattice operators). We find 
the extent to which this hypothesis is satisfied to be surprising, but point 



out in subsec. |4.1| that this was also the case for the AB = 2 operators which 
contribute to B-B mixing, whose matrix elements have been computed by 
several groups. 

The calculation described in this paper is the first evaluation, using lattice 
simulations, of the matrix elements of the operators in eqs.(|^)-(H) between 
i?-meson states. The errors in the results presented in eqs.(P3D-(P^ are 
reasonably small, and until unquenched computations become possible, it 
is not phenomenologically necessary to invest a large effort to reduce the 
statistical errors in this quenched calculation. It is desirable, however, to 
establish, beyond any doubt, that the ground state meson has been isolated 
and to confirm explicitly the validity of the arguments presented in sec. ^ that 
this is indeed the case. For this a similar calculation on a larger statistical 
sample will be required. It would also be very useful to evaluate the one-loop 
contributions to the Wilson coefficient functions in the OPE for inclusive 
decay rates. It is probable that these corrections currently represent the 
largest uncertainty in the inclusive rates. 

A similar study of the matrix elements of the operators eqs.(|35|)-(|38[) 



between A^ states is in progress. This is particularly important in view of 
the discrepancy between the experimental measurement in eq. (|l|) and the 
theoretical prediction in eq. (Q) for the ratio r(Af,)/r(i?rf). This calculation 
will help determine whether the discrepancy is due to baryonic matrix ele- 
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ments being larger than expected from quark models, or whether local quark- 
hadron duality, assumed in the phenomenology, is not valid. Technically the 
baryonic matrix elements of four-quark operators are more complicated to 
evaluate than mesonic ones, since it is not sufficient to generate the set of 
light-quark propagators from an arbitrary lattice point to the origin. The 
results will be presented in a following publication. 

Appendix A - Evaluation of the Dij 

In this appendix we discuss the evaluation of the coefficients {Dy} of eq. (|l^ 
for a generic four-quark operator Oj. This requires the evaluation of one 
loop corrections to the matrix elements of the operators {Oi} in both the 
continuum and lattice schemes: 

A, = eg - C5 , (63) 

where the superscripts C and L refer to the continuum renormalization 
scheme (MS) and the lattice regularization respectively. We evaluate the 
matrix elements between four-quark states at zero momentum, regulating 
the infrared divergences by giving the gluon a small mass A. The coefficients 
{Dij} do not depend on the infra-red regulator, although each of Cg and C^^ 
are separately infrared divergent. 

We distinguish five categories of one-loop corrections to the generic oper- 
ator 

O.^bTfqqtfb, (64) 
where A represents both colour and spinor indices: 

1. Corrections to the external lines. These are proportional to Cp, the 
eigenvalue of the quadratic Casimir operator in the fundamental repre- 
sentation {Cp = 4/3 for the SU{3) colour group), and are independent 
of Fj and Fj. 

2. Vertex corrections to each of the quark bilinears. 

3. Corrections in which the gluon couples to both heavy-quark propaga- 
tors. 
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( catpfforvl 




0^ 


(1) 


-log(AV) + i 


CFbTfqqTfb 


(2) 


-log(A2a2) + l 


b{T^TfT'')q qtfb + bVfq q{T''VfT'')b 


(3, odd) 


21og(A2a2) 


b{T^Tf)qq{tfT^)b 


(4, odd) 


log(A2a2) - 1 


b{TfT^)q g(f f T«)6 + b{T''Tf)q q{T^ff)b 


(5, odd) 


-log(A2a2) + fii 


b{TfT'')qq{T^tf)b 


(5, odd) 


-31og(A2a2) + ^2 


^b{Vfa^"'T'')q q{T''a''^'Vf)b 



Table 3: Values of the contributions to 



4. Corrections in which the gluon couples to one light-quark propagator 
and one heavy-quark propagator (at the other vertex). 

5. Corrections in which the gluon couples to both light-quark propagators. 



In table ^ we present the values of the contributions to the Cf^ in the MS 
scheme. For the operators Oi and O3, i^i = — | and VL2 = 0, whereas for O2 
and O4 fii = 1 and Q2 = 0. For a general choice of Oi, Qi^2 depend on the 
precise definition of the operator b(Tfa^'^)qq{a^^Tf)b and on the choice of 
basis for the 7-matrices in D-dimensions. 

The coefficients {Dij} themselves are presented in table where the Xj's 
are defined as: 

Xi = Ci + cl , (65) 

and the numerical values of the Cj and cf are tabulated in table ^. The 
contributions proportional to the Xj's come from the diagrams in the lattice 
regularization. The components proportional to the Cj would be the results 
if the Wilson formulation of the quark action had been used, and those 
proportional to cf are the additional contributions which result from the use 
of the SW-improved action 0. 

^In the notation of ref . ||l9| , the parameters 5i, and 6s are defined by 

A 1 /•+",4. ^ A4(4-Ai)-(A4-A5) J{l-k^) \ , 
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D 


0^ 


1 (1,2) 


1 


r i 1 A — A 1, 

CpoVfq qVfb 


2 (2) 


1-X2 


b\T VfT°-)qqVfb + bVfqq{T°-VfT°')h 


3 (2) 


— Tq 


b(T'^'y°T^'y°T'')a aVH + bV^a ofr''^°f^^°T'')6 


4 (3, odd) 


— X4 


b{T''Tf)qq{ffT'')b 


5 (4, odd) 


-I-X5 


b{TfT'')q q{VfT'')b + h{T^V^)q q{T''Vf)b 


6(4) 


-X6 


b{Tf^OT'')qq{ff-f°T'')b + b{T^^^Vf)qq{T''^^ff)b 


7 (5, odd) 




b{VfT^)qq{T''Vf)b 


8 (5, odd) 


fi2-a;8 


^b{Vfai"'T'')q q{T''a''^ff)b 


9(5) 




\b{TfYT'')q g(T"7'^f f )6 



Table 4: Values of the contributions to Dj^. 



Appendix B - = 2 operators 

There are many parallels between the calculations of spectator effects in 
inclusive decays, which is the main subject of this paper, and that of the 
matrix elements of the Ai? = 2 operators which contribute, for example, to 
the important process of B^-B^ mixing. Indeed, we have recomputed the 
matrix elements of the AB = 2 operators and compared the results to those 
in ref. as a check on our procedures and programs. The calculation of 
the matching factors are also similar in the two cases, and we have exploited 
this fact as a check on our perturbative calculation. Since we disagree with 
one of the terms in the results of ref. |]19[, we briefly discuss the evaluation 



of the matching factors in this appendix. 

We consider a generic AB = 2 operator of the form bTfqbVfq. The 
evaluation of the various contributions to the -Dj/s for AB = 2 operators 
parallels that of the operator Oi for spectator effects in inclusive decays de- 
fined in eq. (|6^); specifically, as explicitly marked in table H, the coefficients 
corresponding to j = 4, 5, 7 and 8 change sign, whilst the remaining coeffi- 



S =^ /•+"^4^4Ai(A4-A5)(2 + r2Ai)- Al(2 + r2Ai) 



47r2 4AiA2 
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ci = f + e = 17 89 


ci = — 2(1 + m) = —11 53 


C2 = di = 5.46 


= n = 0.73 


C3 = = -7.22 


ci^ = h-2d^ -q = 0.33 


C4 = — c = —4.53 


cl = 


C5 = —di = —5.46 


4 = -n = -0.73 


CQ = d2 = -7.22 


0^ = 4 = 0.33 


C7 = — — 5i, = 4.85 


c;^ = 5, = 0.27 


C8 = 6^ = 2.07 


4 = -y^ + s - 5, + 36l = 10.46 


Cg = 4W7 = -4.84 


= 4w;^ - 2/ - 2m + s + 35/j = -4.88 



Table 5: Values of the coefficients q and cf . Their expressions are given in 
terms of the variables defined in refs. and and two new variables 
6v, Ss- The numerical integrals have been recomputed. Some of them differ 



slightly from the results of |19 . 



cients are the same. Thus from table ^ for spectator effects we deduce that 
the -Dj/s for AB = 2 operators are as given in table ^ The Xj's are defined 
in eq. (^) and tabulated in table ||. 

We are particularly interested in the operator whose matrix elements 
contains the non-perturbative QCD effects for B^-B^ mixing: 

hVfq bf^q = 0l = h'qL hpQL , (66) 

where the label L denotes "left". In this case the different operators Oj of 
table 1^ reduce to 3 independent ones: 

Or = bYqRb'jpQR (67) 
On = 2bqL bqn + 2bqR bqi + IYql b'^pQR + b^QR hpQL , (68) 

as well as Ol itself. Step i) of the matching (see section ^ is now given by 
the relation: 

Of = (1 + + ^DnO^^ + ^DnO^^ (69) 

where 

Dl = -22.4 , Dn = -13.8 and Dr = -3.2 . (70) 
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j (category) 




Of 


1 (1,2) 


2 — Xi 




2 (2) 


I - X2 


b{T''TfT'')qb{Tf)q + h{Vf)qb{T°'VfT'')q 


3(2) 


-X3 


h{T''-f^Tf-i^T'')qh{tf)q + h{Tf)qh{T''-i^tf-i^T'')q 


4 (3, odd) 




h{T^Tf)qh{T''tf)q 


5 (4, odd) 


1 +a:5 


b{TfT^)qb{T''tf)q + b{T^Tf)qb{tfT^)q 


6(4) 


-Xq 


b{Tf^OT'')qb{T^-f^Tf)q + b{T^-f°Tf)qb{Tf^^T^)q 


7 (5, odd) 


—VLi + x^ 


b{TfT^)qb{TfT^)q 


8 (5, odd) 


-VL2 + X8 


^b{Tf(T'"'T'')q b{tfa^'^T'')q 


9(5) 




mTf-f^'T^)qb{tfj^'T'')q 



Table 6: Values of the contributions to the -Djj's for AB = 2 operators. 



For the continuum renormalization scheme we have used MS. In 4 — 2e 
dimensions 

bYcr'^'-qL hpCTf^uqi = (12 - 2e)bYqL hpQL , (71) 

from which we derive that fii + = 5 in this scheme. In considering the 
crossing relations between inclusive decays and mixing, we note that in the 
latter process, one of the b fields destroys a quark (so that 67° = b) and the 
other creates an antiquark (so that 67° = —b). 

The results for and Df^f in eq. ([70|) agree with those in the liter- 
ature |]2B|, whereas that for Dji does not (a similar conclusion was reached 
independently by Gimenez |^). In ref. p6[ the quoted result is = — 5.4Q. 



^We believe that the reason is that in eq. (B.16) of ref. there should be a correction: 



eq. (B.26) should be replaced by 



and that in table 3 



D^j^ = i [s + 4a;^ - 2{l + m)] 



Dj^ = -0.38 for r = 1 
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